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(57) Abstract 



The present invention relates to a method for char- 
acterizing single test samples using techniques generating 
multi dimensional responses from which the components 
of the sample can be identified. The method does not 
require any references and is applicable even on samples 
which are less than the number of components they con- 
tain. 
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METHOD FOR CHARACTERIZING SAMPLES 



The present invention relates to methods for characterizing samples. These methods are i.a. 
used to investigate test samples from a production, patients or samples collected in any 
5 other way. 

Background of the invention 

When a sample is to be characterized for components, the components are generally sepa- 
rated from each other in a first step in order to identified and quantified in a later stage. 
10 However, it is not always possible to separate the components or it may not be motivated 
from a time/cost benefit reason. The samples may then be characterized spectroscopically 
whereby the components are identified by means of their unique spectral responses. 

If one has a collection of samples and is aware of which components they comprise, it is, as 
15 a rule, trivial to determine their concentrations spectroscopically. This is due even if the 
spectral responses of the components overlaps each other. If, however, the components are 
unknown, the problem is muck more complicated . The situation was analysed for the first 
time in detail by the mathematics Lawton and Sylvestre ( Technometrics . 13, 617, (1971)), 
who showed that it is impossible to find an unique solution even for a 2-component system. 
20 In 1990 we developed an experimental method, which partly solved this problem (Kubista, 
Chemometrics and Intelligent Laboratory Systems . 7, 273, (1990)). We then showed that if 
one carried out two spectroscopic measurements on each sample, in stead of one as pre- 
viously used, and the measurements were such that the contribution of the components to 
these measurements had the same distribution of the intensities, but of different magnitude, 
25 then both the spectral responses as well as the concentrations of the components could be 
determined. Mathematically, these measurements are described using the equations: 

r 

A « CV or aj (A) = X c u v i W J = !> 2 — n 

r 

30 B = CDV or bj (A.) = £ c^d.v. (A) j = 1 ,2 ... ji 

wherein A is a matrix comprising spectra of the first type measured on the n samples; B is a 
matrix comprising spectra of the second type measured on the same n samples; C is a 
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matrix comprising the concentrations of the r different components in the n samples; V is a 
matrix comprising the normalized spectra of the components; and D is a diagonal matrix, 
the r diagonal elements of which being the ratios between the responses of the components 
obtained in the two measurements. All spectra are digitalized in m points. We showed that 
5 the concentrations of the components (C), their normalized spectral responses (V) and the 
ratio between their responses obtained in the two measurements (D) could be determined 
only outgoing from the information obtained from the spectra as measured (A and B). We 
further described how the number of components of the samples (r) could be estimated. 

10 One restriction using this method is that the number of components are not allowed to ex- 
ceed the number of samples, which from a practical point of view means that the method 
can not be utilized on smaller series of samples and can not be applied on the whole for 
analysing isolated samples. 



15 Several spectroscopic techniques, such as fluorescence, nmr, etc., can generate 2-dimen- 
sional data described by the equation: 

/(a,/7) = >rX/ / (a)c / / / (/7) 

where the signal, I(a t fi)^ is determined as a function of two variables, a and p, and are the 
20 sum of the contribution of the components in each point, which contribution is proportional 
to their concentrations (Cj) and the products of their (normalized) 1 -dimensional responses, 
I£a) and Out of these responses the components can be identified. In a steady state 
fluorescence spectroscopy I t (a) and 1^/3) are the excitation- and emissions spectra of the 
components and are, as a rule, designated IffA^) and ^ eiw (0Z c J, wherein A ex and A em are the 
25 excitation and emission wavelengths. The shape of an excitation spectra of a pure com- 
pound is, in general independent of the emission wavelength used at the measurement, and 
the corresponding is due for its emission spectrum. The fluorescence signal monitored, if 
necessary after a correction for the inner filter effect (Kubista et al, The Analyst . 1 19, 41 7 
(1994)), is proportional to the concentration of the compound. In a sample containing more 
30 compounds the total signal is the sum of the contribution by each component. As fluore- 
scence is measured in an arbitrary unit, eq. 1 contains a proportionality constant (k). 



The information of the 2-dimensional spectrum I(oc,/3) is insufficient to unambiguously 
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determine the spectral responses of the components. Different approximative ways have 
been suggested but these do not function sufficient satisfactorily even for a 2-component 
mixture (Burdick and Tu, J. Chemometrics . 3, 431, (1989)). 

5 The present invention is a method for analysing isolated test samples, or a couple of test 
samples without using references in such a way that the components can be identified. 



Description of the figures. 

Figure 1 . Emission spectra monitored using different excitation wavelengths using a parallel 
10 polarized light (above, left) and a perpendicularly polarized light (above right), respectively. 
Down to the left the calculated emission spectra of the components are shown, and down to 
the right the calculated excitation spectra of the components are shown. 



Figure 2. A) Excitation spectra registered using different emission wavelengths from two 
15 solutions containing POPOP, dimethyl POPOP, antracene, and diphenyl antracene. B) The 
excitation spectra of the components as calculated. 



Figure 3. A) Emission spectra registered using different excitation wavelengths of two 
solutions containing POPOP, dimethyl POPOP, antracene, and diphenyl antracene. B) The 
20 excitation spectra of the components as calculated. 



Brief description of the invention 

The present invention is a method for analyzing test samples in such a way that its compo- 
nents can be identified without the need for any reference data. The method is based upon 
25 the following four steps: 

1 . The test sample is analyzed using a method generating a 3-dimensional response accord- 

ing to : /(a,Ar) -fj, ("tf Wl (/) > 

wherein r is the number of components contributing to the signal, and I t (a) and and 
30 // y) are the arbitrarily normalized 1 -dimensional responses of the components, which 
responses normally consist of spectral or concentration variations. 

2. The number of components r as the samples contain is estimated. 

3. For each component its 1 -dimensional responses I/a) and I/Jty and I t (y) are determined. 
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4. Out of the responses, the components are identified. 
Detailed description of the present invention 

As the title indicates the present invention relates to a method for characterizing isolated test 
5 samples in a way that makes it possible to identify its components without any need for 
using reference samples. This is done through a strategic design of experiments which 
makes it possible to register a 3-dimensional response being proportional to the concentra- 
tions of the components, and the contribution from each component is the product of its 
specific 1 -dimensional responses: 

Such registration can be carried using certain forms of fluorescence spectroscopy, e.g., by 
means of a time disintegrated monitoring of emission/excitation spectra, i.e., the signal is 
registered as a function of excitation wavelength, emission wavelength, and time: 

15 n*„,*„,t) = i>,/, w„y, w„v, o 

In these cases it is often suitable to gather the concentration of the components ci and the 
time declinations to a time dependent concentration: 

> > o = X <v (OA Vi ) 

20 The time can be time after light pulse (whereby c^t) is proportional to the fluorescence 

declination), time after mixing of e.g., a stop-flow experiment (whereby Cj(t) is the variation 
of the concentration of component i with time), time after treatment, such a photo bleaching 
(selective destruction of certain components using light), chromatographic or other form of 
separation, etc. At the analysis of such data the concentration variation of the components 

25 are calculated, as well as their excitation and emission spectra. It is of interest to note that 
intermediate components which are neither present at the beginning (Cj(0) = 0) or at the end 
( c i(°°) " 0) of the experiment can be identified from its calculated spectra. 

There is a further possibility in varying the polarization of the light: 
30 r 

or, if the phase-modulated light is utilized, the frequency of the modulation: 

I W„ >^v) = t, c t I t (A„ )/,. (A em )/, (v) 
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etc. 



There is further a possibility in varying the outer parameters which influences the concentra- 
tions of the components, such as temperature (pressure, volume, etc.): 

5 '(<*..4..r)-2>,(n/, 

or outer parameters which influence the intensity of the responses of the components, such 
as external magnetic fields (electrical fields, etc.): 

7 ^-^^) = l>^ (^Vi (AO. 

10 <-l 

The spectroscopic technique need not be a fluorescence technique. The method can be car- 
ried out using most techniques which generates 3-dimensional responses, e.g., nuclear 
magnetic resonance spectrometry (NMR) mass spectrometry, etc. It can further be carried 
out using most techniques generating 2 -dimensional responses if the responses of the 
15 components influence external parameters. Finally, the method can be used using a tech- 
nique generating 1 -dimensional responses, as well, but then it is necessary that two external 
parameters are varied simultaneously and that their influence on the responses of the com- 
ponents are independent so that their contribution can be factorized. 

20 The invention requests that at least two data points are determined in each of the 3 dimen- 
sions, i.e.: 

I^a) wherein a u a 2 oc n 1>2 

U$ wherein p„P 2 p m m>2 

Ii(r) wherein Yi,Y 2 Yi n>2 

25 

To determine two data points only in all dimensions are, however, of particular meaning as 
the tolerance of the responses calculated then as a rule is insufficient to be able to identify 
the components. On the contrary it is quite excellent to have two data points only in one of 
the dimensions , e.g, 1=2 (and m»2, and n»2). This exhibits the advantage that the nume- 
30 rical treatment of data is made easier as the responses of the components can be calculated 
using fast algorithms such as Procrustes rotation and GRAM (Kubista, Chemometrics and 
Intelligent Laboratory Systems . 7, 273, (1990); Wilson, Sanches & Kowalski, J. Chemo- 
metrics . 3, 493, (1989)). In the general case when all 1, n, and m are greater than 2, the 
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solution method is much more complicated and thus considerably more time consuming 
(Li wo, et al, Computers Chem .. 21, 89-91, (1997)). Furthermore, it is quite often of interest 
to carry out the experiment in such a way that one of m and n are considerably greater than 
the other. The reason hereto is that it as a rule, is sufficient for the identification of the com- 
5 ponents, to determine one of their 1 -dimensional responses with a high accuracy. 

The invention is not limited to determinations that generates 3-dimensional responses but 
even responses of a higher order can be used. In general it should be satisfying that the 
response is linear and that the contribution from each component shall be the product of its 
10 1 -dimensional responses : 

I{a,^ r> S^) = ^J t Wi (m Wi 

Of course, the higher the dimension is the more time consuming the numerical treatment of 
1 5 the determined data will become. However, with regard to the very fast development within 
the computer area this will hardly be a practical limitation in the future. 

The samples to be analysed shall contain substantially the same components, and these shall 
be present in different, relative concentrations. The samples are analysed in pair using a 2- 
20 dimensional method which provides a response which is proportional to the concentrations 
of the components and the product of the 1 -dimensional responses. This can be expressed 

1=1 

wherein I A (a) and are spectra of the two samples which in the following will be called 
A and B, determined as a function of the variables a and p, r is the total number of compo- 
nents contributing to the spectra, I { (a) and 1/$ are the normalized 1 -dimensional responses 
of the components, and c, A and c, B are their concentrations, respectively. In a steady-state 
30 fluorescence spectroscopy I {a) are the normalized excitation spectra of the components, 
ir(^eJ> an d are their normalized emission spectra, I* m (J> e J: 
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1=1 

The information in these spectra is treated in two steps. First the number of components, r, 
5 is determined, and then the 1 -dimensional responses of the components, 

iraj § ™dira e j- 

When r has been determined, the spectral responses of the components. 
The equations: r 

can be written in matrix form as: 
A=XC A M 
15 B=XC B M 

wherein A and B are matrixes comprising the spectra determined, X is a matrix comprising 
the normalized excitation spectra of the components, M is a matrix comprising their norma- 
lized emission spectra, and C A and C B are diagonal matrixes comprising the concentrations 
of the components. By renormalizing one of X or M, the equation system can be rewritten 
20 as: 

A=XM 
B=XDM 

wherein D is a diagonal matrix comprising the ratios between the concentrations of the 
components (D=C B /C A ). Using A and B, X, M and D can be calculated using known 
25 methods such as Procrustes rotation (Kubista, Chemometrics and Intelligent Laboratory 
Systems . 7, 273, (1990); and GRAM (Wilson, Sanches & Kowalski, J. Chemometrics . 3, 
493, (1989). 

As a summary, the present invention relates to a method for experimentally studying two 
30 samples spectroscopically so that the information present in the experimental spectra is 
sufficient to determine the number of components of the samples (r), their spectral re- 
sponses of the 1 st dimension, I^a), their spectral responses of the 2 nd dimension, I/jfy, and 
the ratios between their concentrations (c, A /C| B ). 
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The most apparent use of the invention is for the analysis of two samples containing 
common components. All components need not be common, but the majority of those con- 
tributing spectroscopically should be in common (Booksh & Kowalski, J. Chemotrics, 8, 
287, (1994)). The number of components is arbitrary and can exceed 2. 

5 

Another use of the invention is to characterize single samples by first dividing them into 
two part samples containing the ingoing components in different proportions. This can be 
accomplished in several ways, e.g., by filtering, extracting, chromatographying dialysing, 
centrifuging, precipitating, splitting the sample by means of an electrical field, etc. Alter- 

10 natively, the original sample can be used as one sample, and an aliquot thereof, which is 
created in such a way that the components are present in other proportions, is used as the 
second sample. This aliquot can be obtained by selectively eliminating certain components, 
e.g., by means of adsorption, precipitation, freezing, distillation, selective decomposing 
(e.g., by light, heat, radio lysis), etc. Another possibility is to create two samples from one, 

15 is to change the conditions for the determination, e.g., by changing the temperature, pres- 
sure, etc. Separation methods, such as different types of chromatography are of interest, as 
the components are separated in space, and one, principally arbitrary number of samples can 
be obtained which can be analysed in pair. Using spectroscopic techniques which generates 
2-dimensional spectra in a fast way, then, furthermore, the detection can be made on-line. 



Another use of the invention is to determine the concentrations of the components in one 
test sample in relation to a standard sample with a high degree of accuracy. The standard 
sample and the test sample are analysed as a pair, and the ratio between the concentrations 
of the components is obtained as the diagonal element of the D matrix. 



2-dimensional spectra wherein one of the dimensions is time, are of particular interest, 
whereby time is related to time after a disturbance such as a relaxation time. Today, there 
are e.g., fluorescence instruments by means of which one can determine complete spectra as 
a function of time after lightening (either directly after lightening using a light pulse, or 
30 indirectly using phase modulation technique). This gives using a as time, and P as wave 
length, the equation system: 



20 



25 



r 
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from which r, I t (t) 9 I/A) and (c^/c®) can be determined. 
5 Example 

The invention will be further illustrated in four examples. 
Example 1 

A sample is characterized using fluorescence spectroscopy, where excitation wave length, 
10 emission wave length, and light polarization are varied (Figure 1). This gives raise to a 3- 
dimensional spectrum according to: 

KA^A^d) = ]>>,./,. (A a y £ (A^I, {a) 

In the example 650 different emission wave lengths (m), 1 1 different excitation wave 
15 lengths (n) and 2 different polarizations (a = Oo, called parallel polarization, and a = 90o, 
called perpendicular polarization) (1), are used. From the response determined, In^A^ A emt 
a), first the number of components (r) is estimated to 2 (using a statistic test and a visual 
inspection of the principal components). Then the component specific responses are cal- 
culated. For this purpose one uses the fact that only two data points were registered in one 
20 of the dimensions (polarization) and rewrote the 3-dimensional response to two 2-dimen- 
sional responses. 

'(^,-L>" = 0°) = ]>>,/,. {A a V, KJ/,. (a = 0°) 



= 90°) = 2>,/ ( . (A„y t (^J/,. (a = 90°) 
25 These can be described using the equation system: 
I°=Xa°M 
I 90 =Xa 90 M 

which can be solved using Procrustes rotation (Kubista, Chemometrics and Intelligent Labo- 
ratory Systems . 7, 273, (1990)). This gave the normalized excitation intensities of the com- 
30 ponents as matrix X{I i (A ex )) > (shown down to the right in Figure 1), the normalized emission 
intensities of the components as matrix M {Ji(A e J) (shown down to the left in Figure 1), and 
the ratios between the responses of components to light of different polarization 
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From the calculated component specific responses, in particular the emission spectra, the 
component could be identified as p-bis[2-(5-phenyloxazolyl)]-benzene (POPOP), and antra- 
cene. Finally, by comparing standard spectra of POPOP and antracene the concentrations 
could be estimated to some micro molars. 

5 

Example 2 

Two solutions containing the dye compounds POPOP, dimethyl POPOP, antracene and 
diphenyl antracene in different proportions were prepared. On these fluorescence excitation 
spectra were monitored at several emission wave lengths. The number of components were 
10 determined to 4 using a statistic test, and the excitation spectra of the components (Figure 
1), emission intensities and the relation between their concentrations in the two samples 
were calculated. 

Example 3 

15 On the same solutions as in Example 1 the fluorescence emission spectra were monitored 
using a number of excitation wave lengths . The number of components was determined to 4 
using a statistic test, and the emission spectra of the components (Figure 2), excitation in- 
tensities and the relation between their concentrations in the two samples were determined. 

20 Example 4 

Characterization of samples containing the dye compound thiazole orange and the polymer 
poly(dG) was made. The samples were analysed in pairs using 2-dimensional fluorescence 
spectroscopy. They contains thiazole orange and poly(dG) in the relation 
[thiazole orange]/[poly(dG)] of 0.05 and 0.025. Neither poly(dG) nor the dye compound is 

25 fluorescentic as such but the fluorescence arises when thiazole orange binds to the polymer. 
The samples were analysed in two different ways. In one analysis, the fluorescence exci- 
tation spectra were monitored at different emission wave lengths. The number of fluorescent 
components were identified to two using statistic tests, and their excitation spectra and 
emission intensities were calculated.. In the second analysis, the fluorescence emission 

30 spectra were monitored using a number of excitation wave lengths. Once again the number 
of components was identified to two, and their emission spectra and excitation intensities 
were calculated. 
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CLAIMS 

1 . A method for characterizing a sample, 
characterized in that 

a) a sample, or pair of samples, is (are) characterized using a monitoring technique such that 
5 a multi dimensional response is generated according to 

I(a,Ar,j....) = £<:,/,. ( a y, (/?)/,. (/)/,- (*).... , 

1=1 

b) the response monitored ir broken down to an orthogonal basset e.g., using a principal 
component division, 

10 c) the number of components (r) in the sample is estimated, 

d) the arbitrary normalized 1 -dimensional responses of the components are calculated. 

2. A method according to claim 1, wherein the number of samples is two and these are ana- 
lysed using a method generating a 2-dimensional response according to 

15 /(*,/?) = £/,.(tf)c ; /, (/?) 
1=1 

and the 1 -dimensional responses of the components and the ratios between their concentra- 
tions in the two samples, (Cj A /Cj B ), is calculated by solving the equation system 

I A (a, ^ = ^1,(^1^/3) 
20 -=> 

l\ a ,^) = ^I i (a)cfl i (jT} 

1-1 

3. A method according to claim 2, wherein the two samples are generated from one sample. 

25 4. A method according to claim 1 or 2, wherein one of the samples is used as a standard 
sample to determine the concentrations of the components in a test sample. 

5. A method according to claim 1, wherein a single sample is amalysed using a technique 
generating 3-dimensional response: 

/(ar, fry) = 2>/ 7 < ^ 7 < W 

and the arbitrary normalized 1 -dimensional responses of the components, %(a) and and 
I t ( y) are calculated. 



30 
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6. A method according to claim 1, wherein a single sample is analysed using a technique 
generating a 2-dimensional response simultaneously as external parameters are varied in 
such a way that the concentration of the components are changed in time: 

5 ' 

and the arbitrary normalized 1 -dimensional responses, I {a) and I/jfy and their changes as to 
concentration c/t) is calculated. 

7. A method according to one or more of claims 1-6, wherein more than two data points are 
10 monitored in only two of the dimensions. 

8. A method according to one or more of claims 1-7, wherein the method generating the 
multi dimensional response is fluorescence or nuclear magnetic resonance method. 

15 9. A method according to one or more of claims 1-8, wherein the variations along, at least 
one of the dimensions, is obtained by varying one external parameter, such as time, elec- 
trical or magnetical field, temperature, modulation, or polarisation 

10. A method according to any of claims 8 or 9, for characterizing a test sample by analys- 
20 ing time dependent emission/excitation spectra, where the time relates to time after excita- 
tion, time after the mixing of thge components, time after a certain treatment of the compo- 
nents, such as chromatographic separation or the similar. 

1 1 . A method according to any of claims 8 and 9 for characterizing a test sample by analys- 
25 ing two time dependencies, in combination with at least some other dependency, such as the 

wave length of the light, where the two time dependencies relates to time after excitation, 
time after the mixing of the components, time after the treatment of the components, such as 
a chromatographic separation. 



30 
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Metod for karakterisering av prover 

Den foreliggande uppfinningen tillhor kategorin metoder for 
karakterisering av prover. Dessa anvands bl.a. for art undersoka 
testprover fran produktion, patienter eller prover som insamlats pa annat 



Uppfinningens bakgrund 

Nar ett prov skall karakteriseras med avseende pa dess komponenter 
separeras i allmanhet forst komponentema fran varandra for art senare 
identifieras och mangdbestammas separat. Det ar dock inte alltid mojligt 
art separera komponentema at, eller sa ar det av tids/kostnadsskal inte 
motiverat. Proverna kan da karakterisera spektroskopiskt, varvid 
komponentema identifieras genom deras unika spektrala responser. 

Om man har en uppsattning prover och vet vilka komponenter de 
innehaller ar det i regel trivialt att bestamma deras koncentrationer 
spektroskopiskt. Detta galler aven om komponentemas spektrala 
responser overlappar. Om komponentema daremot ar okanda ar 
problemet mycket besvarligare. Situationen analyserades forsta gangen i 
detalj av matematikema Lawton och Sylvestre (Technometrics, 13, 1971, 
617), som visade att det ar omojligt att finna en unik losning till och med 
for ett 2-komponentsystem. 1990 utvecklade vi ett experimentellt 
forfarande som delvis loste detta problem (Kubista, Chemometrics and 
Intelligent Laboratory Systems, 7, 273, 1990). Vi visade att om man 
utforde tva spektroskopiska mamingar pa varje prov, istallet for som 
tidigare endast en, och matningama var sadana att komponentemas 
bidrag till dessa hade samma intensitetsfordelning, men olika magnitud, 



satt. 




CP 
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sa kunde bade komponenternas spektrala responser och koncentrationer 
bestammas. Matematiskt beskrivs dessa matningarna med ekvationerna: 



dar Aaren matris som innehaller spektra av det forsta slaget uppmatta pa 
de n st proverna; B ar en matris innehallande spektra av det andra slaget 
uppmatta pa samma n prover; C ar en matris innehallande de r olika 
komponenternas koncentrationer i de n st proverna, V ar en matris 
innehallande komponenternas normaliserade spektra och D ar en 
diagonalmatris vars r diagonalelement ar kvoterna mellan 
komponenternas responser i de tva matningarna. Samtliga spektra ar 
digitaliserade i m punkter. Vi visade att endast utifran informationen i de 
uppmatta spektra (A och B) s£ kunde komponenternas koncentrationer 
(C), deras normaliserade spektrala responser (V) samt kvoten mellan 
deras responser i de tva matningarna (D) bestammas. Vi beskrev ocksa 
hur antalet komponenter i proverna (r) kunde uppskattas. 

En begransningen med detta tillvagagangssatt ar att antalet komponenter 
inte far 6 verstiga antalet prover, vilket i praktiken innebar att metoden 
inte ar tillampbar pa mindre provserier och kan overhuvudtaget inte 
appliceras for analys av enstaka prover. 



A = CV eller a.(i) = X^ v /W 



j=l,2...n 



B = CDV eller bj{A) = ^d^A) 



j=l 5 2...n 



1=1 



Flera spektroskopiska tekniker, sasom fluorescens, nmr etc., kan 
generera 2-dimensionella data som beskrivs av ekvationen: 
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/(ar,/7) = /r£/^)c,. /,(/?) 

1=1 

dar signal en, I(a,JJ) , mats som funktion av tva variabler, a och (3, och ar 
i varje punkt summan av komponentemas bidrag, som ar proportionella 
mot deras koncentrationer (q), och produkten av deras (normaliserade) 1- 
dimentionella responser, /, (a) och /. . Fran dessa responser kan 
komponenterna identifieras. I steady-state fluorescensspektroskopi 
3i I. (a) och L t (jff) komponentemas excitations- och emissionsspektra, 

och betecknas i regel och /; m (^ m ) , dar X cx och X em ar excitations- 

och emissionsvaglangdema. Formen hos ett rent amnes 
excitationsspektrum ar i allmanhet oberoende av den emissionsvaglangd 
som anvands vid matningen, och motsvarande galler for dess 
emissionsspektrum. Den uppmatta fluorescenssignalen, om nodvandigt 
efter korrektion for innerfilter effekten (Kubista et al. ? The Analyst, 1 1 9, 
417, 1994), ar proportionell mot amnets koncentration. For ett prov som 
irmehaller flera amnen ar den totala signalen summan av komponentemas 
bidrag. Eftersom fluorescens mats i godtyckliga enheter, innehaller eq. 1 
en proportionalitetskonstant (k). 

Informationen i det 2-dimensionella spektrumet J r (a, /J) ar otillracklig for 
att entydigt bestamma komponentemas spektrala responser. Olika 
approximativa tillvagagangssatt har fbreslagits, men dessa fungerar inte 
fxillt tillfredsstallande ens for 2-komponentsblandningar (Burdick och 
Tu 5 J. Chemometrics, 3 5 431, 1989). 

Den fbreliggande uppfinningen ar ett fbrfarande att analysera ett enstaka 
testprov, eller par av testprover, utan att anvanda referenser, pa ett sadant 
satt att komponenterna kan identifieras. 



PCT/SE9 8 / 0 140b 

€P 1 4 -08- 1998 



Beskrivning av figurer 

Figur 1. Emissionsspektra uppmatta med olika excitationsvaglangder 
med parallell polariserat (overst till vanster) respekive med vinkelratt 
polariserat (overst till hoger) ljus. Nederst till vanster visas 
komponenternas beraknade emissionsspektra och nederst till hoger visas 
komponentemas beraknade excitationsintensiteter. 

Figur. 2. A) Excitationsspektra uppmatta med olika emissionsvaglangder 
av tva losningar innehallande POPOP, dimetylPOPOP, antracene och 
difenylantracene. B) Komponenternas beraknade excitationsspektra. 

Figur 3. A) Emissionsspektra uppmatta med olika excitationsvaglangder 
av tva losningar innehallande POPOP, dimetylPOPOP, antracene och 
difenylantracene. B) Komponenternas beraknade excitationsspektra. 
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Kortfattad beskrivning av uppfinningen 

Den foreliggande uppfinningen ar ett forfarande att analysera testprover 
pa ett satt att dess komponenter kan identifieras, utan att referensdata 
behovs. Forfarandet bygger pa foljande fyra steg: 

1 . Testprovet analyseras med en metod som genererar en 3-dimentionell 

r ^ 

respons enligt: l{ct,p,y) = YJi ^> C# 7 / (/) > dar r ar antalet 

komponenter som bidrar till signal en och 1. {a) , I. (fi) och I. (/)ar 
komponentemas arbitrart normaliserade 1-dimentionella responser, 
som vanligen utgors av spektrala eller koncentrationsvariationer. 

2. Antalet komponenter, r, som provema innehaller uppskattas. 

3. For vaije komponent bestams dess 1-dimentionella responser /. (a) , 

{JS) och /. {r) . 

4. Fran responserna identifieras komponentema. 

Detaljerad beskrivning av uppfinningen 

Som rubriken antyder ar den foreliggande uppfinningen ett forfarande att 
karakterisera ett enstaka testprov pa ett satt som gor det mOjligt att 
identifiera dess komponentema utan att behova anvanda referensprover. 
Detta sker genom ett strategiskt upplagg av experimented som gor det 
mojligt att registrera en 3-dimentionell respons som ar proportionell mot 
komponentemas koncentrationer, och bidraget fran vaije komponent ar 
produkten av dess specifika 1 -dimentionella responser: 



Sadan matning kan t.ex. utforas med vissa former av 
fluorescensspektroskopi, t.ex. genom tidsupplost matning av 
emissions/excitationsspektra 5 dvs signalen registreras som funktion av 
excitationsvaglangd, emissionsvaglangd och tid: 



r 
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'(4^,0 -!>,/, (^y, (^)/,- (0 

I dessa fall ar det ofta lampligt att sammanfora 
komponentkoncentrationema c { och tidsavklingningarna till en 
tidsberoende koncentration: 

Tiden kan vara tid efter ljuspuls (varvid c,. (/) ar proportionellt mot 
fluorescensavklingningen), tid efter blandning i t.ex. stop-flow 
experiment (varvid c % (t) ar variationen hos komponent i:s koncentration 
med tiden), tid efter behandling, sasom fotobleaching (selektiv 
destruktion av vissa komponenter med ljus), kromatografisk eller annan 
form av separation, etc. Vid analys av dylika data beraknas 
komponentemas koncentrationsvariationer samt deras excitations- och 
emissionsspektra. Det ar intressant att notera att intermediara 
komponenter, som varken ar narvarande i borjan (c, (0) = 0) eller i slutet 
(c,- (oo) = 0 ) av experimented kan identifieras utifran dess beraknade 
spektra. 

Det ar ocksa mqjligt att variera ljusets polarisation: 

= ^1, i^Vi W«)/ f (*) 
eller, om fasmodulerat ljus anvands, modulationens frekvens: 

/<^^.*) = l>,'i i^exVi <4J/, (") 
osv. 

Det ar ocksa mqjligt att variera yttre parametrar som paverkar 
komponenterna koncentrationer, sasom temperatur (tryck, volym, etc): 



rUI / bt y b / U '! 4 o b 
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r 



/(A 



eller yttre parametrar som paverkar intensiteten i komponentemas 
responser, sasom externa magnetfalt (elektriska fait, etc) 



Den spektroskopiska tekniken behover naturligtvis inte vara fluorescens. 
Forfarandet kan utforas med de fiesta tekniker som genererar 3- 
dimensionella responser, t.ex. karnmagnetisk resonans (nmr), 
masspektrometri etc. Den kan aven utforas med de fiesta tekniker som 
genererar en 2-dimensionell responser om komponentemas responser 
paverkas externa parametrar. Slutligen kan forfarandet aven anvandas 
med en teknik som genererar 1-dimensionell respons, men da kravs att 
man samtidigt varierar tva olika externa parametrar och att dessas 
inverkan pa komponentemas responser ar oberoende sa att deras bidrag 
kan faktoriseras. 

Uppfinningen kraver att atminstone tva data punkter bestams i var och en 
av de 3 dimensionerna, dvs: 
L t (a) dar a l9 a 2 ..xz m 1>2 

/,(/?) dar^.A-A m>2 

dar y x ,r 2 n>2 
Att enbart bestamma tva datapunkter i samtliga dimensioner ar dock 
sallan meningsfullt, eftersom noggrannheten i de beraknade responsema 
ar da i regel otillracklig for att komponentema ska kunna identifieras. 
Daremot gar det utmarkt att endast ha tva datapunkter i en av 
dimensionerna, t.ex. 1=2 (och m»2 och n»2). Det har dessutom 



r 



7(4 



t-M) = ZCJ; (M) . 
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fordelen att den numeriska behandlingen av data underlattas, eftersom 
komponentemas responser kan beraknas med snabba algoritmer sasom 
Procrustes rotation och GRAM (Kubista, Chemometrics and Intelligent 
Laboratory Systems, 7, 273, 1990; Wilson, Sanches & Kowalski, J. 
Chemometrics, 3 493, 1989). I det generella fallet nar samtliga 1, n och m 
ar storre an 2, ar losningsforfarandet mer komplicerat och darmed 
betydligt mer tidskravande (Liwo et al., Computers Chem. 89-96, 21, 
1997). Dessutom, ar det ofta aktuellt att utfora experiment sa att endera 
utav m och n ar betydligt storre an den andra. Skalet ar att det i regel ar 
tillrackligt for komponentemas identifiering, att med hog noggrannhet 
bestamma en utav deras 1-dimensionella responser. Detta ar t.ex. fallet i 
exempel 1 , dar endast komponentemas emissionsresponser bestamds 
med hog noggrannhet. 

Uppfmningen ar inte begransad till matningar som genererar 3- 
dimensionell respons, utan aven hogre ordningens responser kan 
anvandas. Generellt galler att responsen skall vara linjar och bidraget 
fran varje komponent skall vara produkten av dess 1-dimensionella 
responser: 

I{a,/3, r ,S....) = {ay, (J3)I t (r)/, (S).... 

Naturligtvis, ju hogre dimensionen ar ju mer tidskravande blir den 
numeriska behandlingen av matdata. Dock, med tanke pa den mycket 
snabba utveckling som sker inom dataomradet, kommer detta knappast 
vara en praktisk begransning i framtiden. 



Provema som analyseras skall innehalla mestadels samma komponenter, 
och dessa ska forekomma i olika relativa koncentrationer. Proverna 
analyseras parvis med en 2-dimentionell metod som ger en respons som 
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ar proportionell mot komponentemas koncentrationer och produkten av 
deras 1 -dimentionella responser. Detta kan uttryckas: 

1=1 

dar I A {a,P) och I B {a,j3) ar spektra pa de tva provema, som i 
fortsattningen betecknas A och B, uppmatta som flinktion av variablema 
a och p\ r ar det totala antalet komponenter som bidrar till spektra, /, (a) 
och l,(jff) ar komponentemas normaliserade 1 -dimentionella responser, 
och c, A och cf ar deras koncentrationer. I steady-state 
fluorescensspektroskopi ar /, (a) komponentemas normaliserade 
excitationsspektra, och /, {fi) ar deras normaliserade 

emissionsspektra, /,"(>*«.) : 

Informationen i dessa spektra behandlas i tva steg. Forst bestams antalet 
komponenter, r, och sedan komponentemas 1 -dimensionella responser, 

I?(A a ) och irw„). 

Nar r bestamds beraknas komponentemas spektrala responser. 
Ekvationema: 

1=1 
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kan i matrisform skrivas: 

A=XC A M 
B=XC B M 

Dar A och B ar matriser innehallande de uppmatta spektra, X ar en matris 
innehallande komponenternas normaliserade excitationsspektra, M ar en 
matris innehallande deras normaliserade emissionsspektra, och C A och 
C B ar diagonalmatriser innehallande komponenternas koncentrationer. 
Genom art omnormalisera endera X eller M, kan ekvationssystemet 
omskrivas till: 

A=XM 
B=XDM 

Dar Daren diagonalmatris innehallande kvoterna mellan 
komponenternas koncentrationer (D=C B /C A ). Utifran A och B kan X, M 
samt D beraknas med kanda metoder sasom t.ex. Procrustes rotation 
(Kubista, Chemometrics and Intelligent Laboratory Systems, 7, 273, 
1990) och GRAM (Wilson, Sanches & Kowalski, J. Chemometrics, 3 
493, 1989). 

Sammanfattningsvis utgors den foreliggande uppfmningen av ett 
forfarande att experimentellt studera tva prover spektroskopiskt sa att 
den information som finns i de experimentella spektra ar tillracklig fbr 
att bestamma antalet komponenter i proverna (r), deras spektrala 
responser i den 1 :a dimensionen, /, (a) , deras spektrala responser i den 
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2:a dimensionen, /,. (/?) , samt kvotema mellan deras koncentrationer 
(c, B /c, A ). 

Den mest uppenbara anvandning av uppfinningen ar for analys av tva 
prover som innehaller gemensamma komponenter. Alia komponenter 
behover inte vara gemensamma, men majoriteten av dem som bidrar 
spektroskopiskt bor vara det (Booksh & Kowalski, J. Chemotrics, 8, 287, 
1994). Komponentemas antal ar godtyckligt, och kan overstiga 2. 

En annan anvandning av uppfinningen ar att karakterisera enstaka prov 
genom aft forst dela upp det i tva delprover som innehaller de ingaende 
komponenterna i olika proportioner. Detta kan astadkommas pa ett flertal 
satt, t.ex. genom att provet filtreras, extraheras, kromatograferas, 
dialyseras, centrifugeras, utfalls, delas upp med elektriskt fait, etc. 
Alternativt kan ursprungsprovet anvandas som ett prov och en delmangd 
av detta, som skapas sa att komponenterna fbrekommer i andra 
proportioner, anvands som det andra provet. Denna delmangd kan 
astadkommas genom att selektivt avlagsna vissa komponenter t.ex. 
medelst adsorption, utfallning, nedfrysning, destination, selektiv 
degradering (t.ex. med ljus, varme, radiolys), etc. Ytterligare en 
mqjlighet att skapa tva prover fran ett ar att andra betingelsema for 
matningen, t.ex. genom att andra temperatur, tryck etc. 
Separationsmetoder, sasom olika former av kromatografi, ar intressanta, 
eftersom komponenterna separeras rumsligt och ett, i princip, godtyckligt 
antal prover kan erhallas som kan analyseras parvis. Med 
spektroskopiska tekniker som snabbt genererar 2-dimensionella spektra 
kan dessutom detektionen ske on-line. 
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Ytterligare en anvandning av uppfinningen ar att med hog noggrannhet 
bestamma koncentrationerna av komponenterna i ett testprov relativt ett 
standardprov. Standardprovet och testprovet analyseras som par, och 
kvoten mellan komponentemas koncentrationer erhalls som 
diagonalelementen i D matrisen. 

Av sarskilt intresse ar 2-dimensionella spektra dar en av dimensionerna 
ar tid, avseende tid efter storning, sasom relaxationstid. I dag firms, Lex., 
fluorescensinstrument med vilka man kan mata kompletta spektra som 
funktion av tid efter belysning (antingen direkt efter belysning med 
ljuspuls, eller indirekt med fasmodulationsteknik). Detta ger med a som 
tid och (3 som vaglangd ekvationssystemet: 



fran vilket r, /. (/) , /. (A) och (cf/cf ) kan bestammas. 

Exempel 

Uppfinningen illustreras med fyra exempel 
Exempel 1 

Ett prov karakteriseras med fluorescensspektroskopi dar 
excitationsvaglangd, emissionsvaglangd och ljusets polarisation varieras 
(Figur 1). Detta ger upphov till ett 3-dimensionellt spektrum enligt: 



r 



I A (t,A) = Y J I i U)cfI i (A) 



I B (t,A) = f j I i (t)cfl i (A) 



r 



7(4 



^.«)=E*,/i(*.yi(4.)/,<«) 
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I exemplet anvands 650 olika emissions vaglangder (m), 1 1 olika 
excitationsvaglangder (n) och 2 olika polarisationer (a = 0° benamnt 
parallell polarisation och a = 90° benamnt vinkelratt polarisation, 1). 
Fran den uppmatta responsen, I a), uppskattades forst antalet 

komponenter (r) till 2 (med statiskt test och visuell inspektion av 
principalkomponenterna). Sedan beraknades de komponentspecifika 
responserna. For detta utnyttjade man det faktum att endast tva 
datapunkter registrerats i en av dimensionerna (polarisation), och skrev 
om den 3-dimensionella responsen till 2 st 2-dimensionella responser: 

= 0°) = X^A (^Vi W-Vi ( a = °°) 
1=1 

M„ = 90°) = £ c,J, U„ V, U em )/, {a = 90°) 

1 = 1 

Dessa kan beskrivas med ekvationssystemet: 

I°=Xa°M 

I 90 =Xa 90 M 

som kan losas med Procrustes rotation (Kubista, Chemometrics and 
Intelligent Laboratory Systems 7, 273, 1990). Detta gav komponenternas 
normaliserade excitationsintensiteter som matris X(7 i visas nederst 
till hoger i figur 1), komponenternas normaliserade emissionsintensiteter 
som matris M (I. {Z em ) , visas nederst till vanster i figur 1), samt kvotema 
mellan komponenterna responser av ljus med olika polarisation 
(/, {a = 90°)//, {a = 0°)) = 0.25 och (/, {a = 90°)/7 2 (a = 0°))= 0.94). Fran de 
beraknade komponentspecifika responserna, sarskilt emissionsspek-tra, 
kunde komponenterna identifieras som p-bis[2-(5-phenyloxazolyl)] 
benzen (POPOP) och antracen. Slutligen , genom att jamfbra med 
standardspektra av POPOP och antracen, kunde koncentrationerna 
uppskattas till nagon mikromolar. 
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Exempel 2 

Tva losningar innehallande fargamnena POPOP 3 dimetylPOPOP 5 
antracene och difenylantracene i olika proportioner tillreddes. Pa dessa 
mattes fluorescensexcitationsspektra vid ett flertal emissionsvaglangder. 
Komponenternas antal bestamdes till 4 med statiska test, och kompo- 
nenternas excitationsspektra (Figur 1 ), emissionsintensiteter samt for- 
hallande mellan deras koncentrationer i de tva proverna beraknades. 

Exempel 3 

Pa samma losningar som i exempel 1 mattes fluorescens emissions- 
spektra vid ett flertal excitationsvaglangder. Komponenternas antal 
bestamdes till 4 med statiska test, och komponenternas emissions- 
spektra (Figur 2), excitationsintensiteter, samt forhallande mellan deras 
koncentrationer i de tva proverna beraknades. 

Exempel 4 

Karakterisering av prover innehallande fargamnet tiazolorange och 
polymeren poly(dG). Proverna analyseras parvis med 2-dimensionell 
fluorescensspektroskopi . De innehaller tiazolorange och poly(dG) i 
forhallandena [tiazolorange]/[poly(dG)] 0.05 och 0.025. Varken 
poly(dG) eller fargamnet ar fluorescenta i sig sjalva, utan fluorescensen 
uppkommer nar tiazolorange binder till polymeren. Proverna analysera- 
des pa tva oberoende satt. I en analys mattes fluorescensexcitations- 
spektra med olika emissionsvaglangder. Antalet fluorescenta compo- 
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nenter identifierades till tva med statiska tester, och deras excitations- 
spektra och emissionsintensiteter beraknades. I den andra analysen 
mattes fluorescensemissionsspektra med ett antal olika excitations- 
vaglangder. Ater identifierades antalet komponenter till tva, och deras 
emissionsspektra och excitationsintensiteter beraknades. 
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Patentkrav 

1. Ett forfarande att karakterisera testprover som kannetecknas darav att: 

a) ett prov, eller par av prover, karakteriseras med matteknik sa att 
en multidimensionell respons genereras enligt: 

I{a,J3, r ,S....) = £>,./,. W 0*> 7 < " ' 

b) den uppmatta responsen bryts ned till ett ortogonalt basset, t.ex. 
med principalkomponentuppdelning, 

c) antalet komponenter i provet (r) uppskattas, 

d) komponenternas arbitrart normaliserade 1-dimensionella, 
responser beraknas. 

2. Ett forfarande enligt krav 1 dar proverna ar tva till antalet och 
analyseras med en metod som genererar 2-dimentionell respons 

enligt: I(a y /7) = £ {a)c i I i {/J) och komponenternas 1-dimensionella 
i=i 

responser och kvoterna mellan deras koncentrationer i de tva 
proverna, c, B /cf , beraknas genom att losa ekvationssystemet 

3. Ett forfarande enligt krav 2 dar de tva proverna genereras utifran ett 
prov. 

4. Ett forfarande enligt krav 1 eller 2 dar ett av proverna anvands som 
standardprov for att bestamma komponenternas koncentrationer i ett 
testprov. 

5. Ett forfarande enligt krav 1 dar ett enstaka prov analyseras med en 
teknik som genererar 3-dimensionell respons: 

/(«, Ar)- (artf 
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och komponenternas arbitrart normaliserade 1 -dimensionella 
responser, I. (a), 7 i {JJ) oz\il t (/)beraknas. 
6. Ett forfarande enligt krav 1 dar ett enstaka prov analyseras med en 
teknik som genererar 2-dimensionell respons samtidigt som externa 
parametrar varieras sa att komponenternas koncentrationer andras i 
tiden: 



och komponenternas arbitrart normaliserade 1 -dimensionella 
responser, 7. (a) , 7 i {/J) och deras koncentrationsforandringar c t (t) 

beraknas. 

7. Ett forfarande enligt nagot av kraven 1 till 6 dar mer an tva 
datapunkter registreras i endast tva av dimensionerna. 

8. Ett forfarande enligt nagon av kraven 1 till 7 dar metoden som 
genererar den multidimensionell responsen ar fluorescens eller 
karnmagnetisk resonans. 

9. Ett forfarande enligt nagot av kraven 1 till 8 dar variationerna utmed, 
atminstone en av dimensionerna, astadkoms genom att variera en 
extern parameter, sasom tid, elektriskt eller magnetiskt fait, 
temperatur, modulation, polarisation 

10. Ett forfarande enligt nagot av kraven 8 och 9 att karakterisera ett 
testprov genom att analysera tidsberoende 

emissions/excitationsspektra, dar tiden avser tid efter excitation, tid 
efter komponenternas blandande, tid som komponenterna utsatts for 
viss behandling sasom kromatografisk separation och liknande. 

1 1 . Ett forfarande enligt nagot av kraven 8 och 9 att karakterisera ett 
testprov genom att analysera tva tidsberoenden, i kombination med 
atminstone nagot annat beroende, sasom ljusets vaglangd, dar de tva 



r 



/(*,A/) = £c f (/)/,. (<)/,. 00) 
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tidsberoenden avser tva utav tid efter excitation, tid efter 
komponentemas blandande, tid som komponenterna utsatts for 
behandling sasom kromatografisk separation 
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Sammanfattning 

Uppfinningen ar ett forfarande att karakterisera enstaka testprover med 
tekniker som genererar flerdimensionella responser fran vilka provets 
komponenter kan identifieras. Forfarandet kraver inga referenser och ar 
tillampbart aven pa prover som ar farre till antalet an det antal 
komponenter som de innehaller. 
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V. Resoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 



1 . Statement 



Novelty (N) Claims 1-11 YES 

Claims NO 

Inventive step (IS) Claims 1-11 YES 

Claims NO 

Industrial applicability (IA) Claims 1-11 YES 

Claims NO 



2. Citations and explanations 

The invention relates to a method for spectroscopic (e.g. 
fluorescence or nuclear magnetic resonance) characterization 
of a single sample using a monitoring technique that generates 
a multi-dimensional response proportional to the concentration 
of different components in the sample. The response monitored 
is broken down to an orthogonal base, whereafter the number of 
components in the sample is estimated and arbitrary normalized 
1-dimensional responses of the components are calculated. 

The invention also includes the idea of varying an external 
parameter, such as time, electric or magnetic field, 
temperature, • modulation or polarisation in order to accomplish 
variations along at least one of the dimensions. 

In prior art methods, which could be represented by the cited 
article "A New Method for the Analysis of Correlated Data 
Using Procrustes Rotation which is Suitable for Spectral 
Analysis", Chemometrics and Intelligent Laboratory Systems, 7 
(1990), the number of components are not allowed to exceed the 
number of samples. From a practical point of view this method 
can not be utilized on small sample series and can not be 
applied to the analyzation of isolated samples. 

The claimed method does not require any references and is 
applicable even on samples that are fewer than the number of 
components they contain and also when the number of components 
are unknown.^ 

Therefore, the invention according to claims 1-11 is novel, is 
considered to involve an inventive step and to be industrially 
applicable. 
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